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Abstract 

We derive thermodynamic Bethe ansatz equations describing the vacuum energy of the 
SU (2N) I Sp(N) nonlinear sigma model on a cylinder geometry. The starting points are the 
recently-proposed amplitudes for the scattering among the physical, massive excitations of 
the theory. The analysis fully confirms the correctness of the S-matrix. We also derive closed 
sets of functional relations for the pseudoenergies (Y-systems). These relations are shown to 
be the k — > oo limit of the Sp(k + l)-related systems studied some years ago by Kuniba and 
Nakanishi in the framework of lattice models. 
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1 Introduction 



Since the early days of modern particle physics, the interest in four dimensional non-Abelian 
gauge theories has also strongly encouraged the study of two dimensional nonlinear sigma mod- 
els. The main, well-known reason is that the two theories share the remarkable properties of 
renormalizability, confinement and asymptotic freedom. Besides the obvious advantage of work- 
ing in a lower dimensional setup, many 2D nonlinear sigma models exhibit the even more useful 
property of integrability. Exact solubility offers the unique opportunity to investigate non pertur- 
bative phenomena in quantum field theory and statistical mechanics via non perturbative tools 
like the exact S-matrix [1], the Form Factor [2] and the thermodynamic Bethe ansatz [3] (TBA) 
methods. Along the lines of the earlier work by Fateev and Zamolodchikov [4] but also motivated 
by new, interesting applications in condensed matter physics and string theory [5] there has been 
recently a renewed interest in 2D integrable sigma models [6-13] . 

When dealing with an integrable quantum field theory, a well established strategy is to derive 
the scattering matrix by solving suitably imposed Yang-Baxter equations, together with the 
unitary and crossing constraints. The S-matrix, via the so-called bootstrap procedure, gives 
access to the full infrared mass spectrum of the theory and via the Form Factors, to the large 
distance behaviour of correlation functions. 

In order to explore the ultraviolet phase of the model and, at the same time, to check the 
correctness of the S-matrix, the standard tool is the thermodynamic Bethe ansatz. The TBA 
method is susceptible to all the details of the scattering matrix, it uses the full mass spectrum and 
it is sensitive to the ambiguous (CDD) scalar factors contributing to the amplitudes. It leads to 
numerically solvable integral equations for the vacuum energy on a cylinder geometry and exact 
values for central charge and scaling dimensions of the conformal field theory (CFT) describing 
the model in its deep ultraviolet limit. 

Following the idea of Zamolodchikov [14], and starting from a given TBA equation, it is 
possible to derive sets of functional relations called Y-systems. The Y-systems are "universal" in 
the sense that there is a one-to-many correspondence with TBAs. A Y-system can be transformed 
by imposing different analytic properties, to physically distinct TBAs. (The relationship between 
Y-systems and TBA equations is practically the same as the one between differential equations 
and their solutions.) 

Many examples of functional relations which exhibit elegant Lie algebraic-related structures 
are now known. Most of them are, or were at some stage, purely conjectures (see, for example, 
[15-25]). An important family of such systems is encoded on a pair of Dynkin diagrams of ADE 
type [20-22]; a second family is in correspondence with all the simple Lie algebras 4 and it was 
studied in the works by Kuniba and Nakanishi [17] in the closely related framework of integrable 
lattice models. Many, but by no means all, solutions of these systems have been successfully 
related to perturbations of minimal CFTs and to sigma models [3,4,6,8,16-27]. 

In this paper, such an example of TBA system is derived in the context of integrable non- 
linear sigma models. Being infinite in one direction, the equations require a truncation for the 
ultraviolet analysis. The truncated systems are found to describe the perturbation of the coset 
SU (27V)fc + i / 'Sp(N)k+i via the relevant operator (l,adj). 

The starting point is the recently proposed [12] S-matrix for the sigma models defined on the 
symmetric space AIL (This space is isomorphic to the group coset manifold SU (2N) / Sp(N) .) 
The S-matrix is of SU(2N) Gross-Neveu type, with CDD factors modified in such a way that their 
poles describe only the particle multiplets related to the even fundamental weights of SU(2N). 

4 These two families of Y-systems are partially overlapping. 
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The plan of the paper is as follows. The model is briefly introduced in section |2 The 
derivation of the TBA equations is sketched in section 01 The Y-system is derived in section 0] 
and its equivalence to a particular limit of the C^+i-related 5 system of [17] is discussed in section|21 
Using the dilogarithm identities of [17] the ultraviolet central charge is found to be c = 2N 2 —N—1: 
a result which is fully consistent with the interpretation of the model as a perturbation of the 
SU(2N)/ Sp(N) conformal field theory. Section El contains our conclusions, while the explicit 
definitions of the kernel functions are collected in the appendix at the end of the paper. 



2 The All sigma models 

The classical action of the quantum field theory in question is defined through the functional 

%] = ^ j d 2 xTr[d^g{x)d^ l {x)} , (1) 

where g(x) belongs to a matrix representation of SU(2N) constrained on the symmetric space 
SU(2N)/ Sp(N). Then the action Snlsm[4>] = S[g(4>)] of the model is obtained from S[g] by 
solving the constraint in terms of the independent set of fields {</>} and substituting the resulting 
g((p{x)) into ©. 

The spectrum of the All sigma models consists in multiplets transforming in the represen- 
tations with even fundamental highest weights // 2 a of SU(2N). The masses of these multiplets 
are 

m 2a = Msm(^), a=l,...,N-l, (2) 



\N 

(M = M(A) is the mass scale). The amplitude 26 f° r scattering of a particle 2a on the particle 
2b is obtained by fusion of #2,2- The latter was proposed in [12], and has the form [28] 

S2a,2b{9) = X 2a ,2b(8) cr 2a,2b{9)S2a,2b(9) • (3) 

S2a,2b(@) is the SU(2N) invariant R-matrix spectral decomposed into its irreducible represen- 
tations [29], (T2a,2b(Q) is the cross-unitarity factor and X 2ai2 &(#) is again obtained by fusion of 
X 2)2 (0). They have the following forms [28,29] 

X2a,2b(0) = f[Y[X{6 + 4iriA(i + j-l-?±±)) 

i=l J=l 

= ex P {-/_^e^(ir o >)-^)} , (4) 

with g = 2N, A = ±, and 

X ^ = - Ufa A 'A /O ' °WbiP) = eX P i / — e J ^2a,2b (w) } • 5) 

smh((0 - 4vrzA)/2) [7_ oc u J 

The explicit expressions of A^^ftj) an d i ^2a,2b( UJ ) are given in the appendix at the end of the 
paper. (Notice that the factor X2 a ,2b{@) is the same as X a b(9) in the SU(2N) invariant Gross- 
Neveu model of [29], however the A^JVl(oj) are different. This is because an N-dependent form 
To avoid any confusion: by Sp(n) we mean the compact real form for the classical Lie algebra C n . 
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of the Fourier transform is used in this paper.) Using these expressions we find 

<x 2a , 2b (#)X 2a , 2fe (0) = exp | - y°° ^-e^ (y 2a , 2fe M - <U) } , (6) 

^ e^l sinh(2|a;| min(a, b)) sinh(2a;(A r — max(o, b))) 

coshw smh(2a;iV) 

3 Diagonalisation of the transfer matrix and the TBA equations 

The standard derivation of the TBA equations starts by imposing periodicity conditions on 
the wavefunction describing a gas of particles with rapidity set {9i} and confined on a ring of 
circumference C. The quantisation condition is 

<>"-. c ^"' T** i (0.\e j+1 ,...,e // ,e 1 ,...,e j - 1 ) = i , (8) 

T 2a i is the transfer matrix and it is constructed out of the S-matrix as 

T 2a > %, 0i, 9j-i) = -Tr [5 2 „,,2a 1 {Oj ~ &l) ■ ■ ■ S^ZaM ~ M] , (9) 

where the trace is over the representation of the particle of type 2a j. The main problem is then 
to find the largest eigenvalue T$ of T. The solution is [30] 

N M2aj (a _ ( 2a ^ _l • a\ 

T? j = n x^vi-Oi) -OA* n - v ia) ) ■ ( io ) 

i=1)¥j a=1 \0j-Va 3 -ITtAJ 

In equations the set of parameters {i>i } is solution of the Bethe equation [29,30] 



11 I -.(6) /) ,^A,, „ / 11 11 ..(6) » „^A~ „ I " 1 j 



j=l ~ u 3 ' <'^ L ^ JJ a,j "6/ c =l ,3=1 \ "a ~~ "fl 

w& and at, are, respectively, the fundamental weights and roots of SU(2N), M c is the number of 
solutions of type c, and Q,a is a phase factor which is irrelevant for the current calculations. In 
the thermodynamic limit 

M -f oo , £ -» oo , AT/£ = finite , (12) 



the solutions of (jllj) group into complex "strings" of length p and real part v, 



(6) 

op 



=ug+i7rA(p+l-2A;), fc = l,2,...,p. (13) 

Using (|11() in equation (|SJ) with T replaced by To, introducing "magnon" densities Pp(v) for the 

real parts of Va \ densities cr 2a (0) for the particles, and after writing the log of equations ((SJ and 
(fTTj) in terms of the densities we obtain 

N-l oo 

c?2a(e) + a 2a (6) = ^ coshtf + £ F 2a)26 * <T 2b (0) - ]T a p * pf{9) , (a = 1, . . . , N - 1) (14) 

6=1 p=l 
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oo 2JV-1 

a p *a c (e)=P C p + Y, E i M) ^r*# - (c = l,...,2iV-l) . (15) 

9=1 6=1 

(The symbol * denotes the standard convolution.) In 151 ) a 2a and complement the densities 
of particles to the total densities of states, and are called densities of holes. The expressions for 
the kernels Ap^\ K^P and a p , together with some useful identities, are written in the appendix. 
Solving the latter equation with respect to p h q and plugging the solution into equation l|14|): 

JV-l oo 

ft{0) = ^ooahfl-^E^ * 43 (a = l,...,N-l) (16) 

6=1 9=0 

2JV-1 oo 

w = - EE ^ * 4f } * *5w . ( C = i, . . . , 2iv - 1) (I?) 

6=1 g=0 



where p = 1, 2, . . 
and 



c c — c ~c c — c — C C / 1 o\ 

s o = a o > s o = °"o ) s p = Pp i s p = Pp > (4°J 



2i+l _ ~2i+l 



ofg* +i = . (19) 



The new kernel has the following Fourier representation 



z > rV )^V)^ 1 - i - #R j. (20) 

From equation (|17|) we see that the excitations 3? a with p = are massive, whereas those with 
p > 1 are of magnonic type. The next step in the derivation is also standard: one looks for the 
minimum of the free energy F = E — TS subject to the constraint (|17|l.( T is the temperature, E 
and S are, respectively, the energy and the entropy and are expressible in terms of the densities 
s,s.) The final equations are 

N-l oo 

Ll\{6) = Rm 2a coshfl + £ £ ^ * Z ( Q f * Lf_(6) , (a = 1, . . . , N - 1) (21) 

6=1 9=0 

JV-l oo 

^+W = EE^2 iV) *4" ) *4-W. (c=l,...,2iV-l) (22) 

6=1 9=0 

where L£ ± = In (l + e ±e ?), i? = 1/T and e a p = 3£/s£ are the dressed energies (pseudoenergies). 
The free energy per unit length f(R,M) is 

1 N-l r-OO 

f(R,M) = -—J2 d9 m 2a cosh 9Ll a _(e) . (23) 

In the ultraviolet limit, i?M — > 0, /(-R, M) can be related to the central charge of the correspond- 
ing conformal field theory by 

f(R,M)~--^c uv ■ (24) 
For the All sigma model, the natural expectation is 

c uv = Dim[SU(2N)) - Dim[5p(iV)] = 2N 2 - N - 1 . (25) 
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4 The Y-system 

A Y-system is a set of functional relations among the unknown functions 

Y p a {9) = e £ i (e) . (26) 

To derive a Y-system from the TBA obtained in the previous section one starts from eq. ()21|) in 
Fourier representation, multiply it by K^ N \2uo) and use properties (|48|) and (|49|) to get 

£ x20 (2w) £» („) = £ ^^(o,) (2cosh(o,)^ 2a + LJ-H , (27) 

cos \u) 6=1 c=1 q=0 v 



with 



ff/gV) = (l 5 ^ r -\ ( 8 q0 - ) . (28) 

0q { ; V 2 cosh 2 (w); V 2cosh(a;) ' 



Equation (|27|l can also be written as 



AT-l 



cosh(2o;) ?2a 1 t( n -1) f2b 1 \ 

^sh(^ L0+(W) " 2 C Osh(^ ^ J -6 *<«H-M 



6=1 

2 AT-l 



^T^ a M TTT F C N ~ 1] L\ («) ~ Lf ( W ) . (29) 

cosh(w) v ; 2cosh(u;) ^ a6 1 



6=1 

By using L£ + - = e£: 

AT-l 2A-1 

2cosh(2u,)g 2 Q » + 2cosh( W )L 2 « ( w ) = £ i£°^.(a;) - £ i£" > (30) 

6=1 6=1 

or, equivalently 

Y o 2a (# + 2nrA)Y o 2a (0-2^A) = ( 1 + -— - , - ) (l ' 



Y^ 2 " (0 + i7rA) / V ^i 2 " (8 ~ . 

iV-l T (N-i) 2N-1 , v _ 7 (2iv-i) 

n(i+nrw)" n (1+^) -< 31 ) 



(a = 1, . . . , N— 1.) The same procedure applied to equation (|22j) leads to a second set of functional 
relations of SU(2N) Gross-Neveu type [8,16]: 

2N-1 l{ 2N-l) 00 f . 

Y p %9 + ivA)Y p %9 - ivrA) = J] (l + *?(*)) II 1 + > ( 32 ) 

6=1 <?=0 V 9 

where p = 1, 2, . . . , 00, c = 1, 2, . . . , 2iV — 1 and by definition Y~ 2l+1 = 00. The Y-system (j3U32j) is 
pictorially represented in figure ^ There is a one-to-one correspondence between the Y functions 
and the nodes in the diagram. Black nodes correspond to massive pseudoenergies, while white 
ones are of magnonic type; a link connecting a node (p, a) and a node (p', a') indicates that when 
Y p a appears on the l.h.s. of equation (or there is a term containing YS on the r.h.s. 
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• • • • 



Figure 1: The TBA diagram for the SU(6)/ Sp(3) sigma model. 



The kind of term on the r.h.s. is different for different styles, horizontal and vertical links. Arrows 
indicate one-way directed links. 

The ultraviolet central charge can then be written in terms of the sets of numbers {3^} 
and ^-independent solutions of the full Y-system and of the subsystem obtained after 

eliminating all the massive nodes: 



6 

cuv = — n 



E 



all nodes 

where L is the Rogers dilogarithm function. 



E 

magnonic nodes 



L 



1 



(33) 



5 Truncation of the Y-system and its ultraviolet analysis 

One of the main results of [17] is an exact conjecture for sum-rules of the type (|3*3*|) for Y-systems 
related to the classical Lie algebras 



- y 



L 



a[KN] 



7T 



all nodes 



i + y a P [KN] 



l Dim [G(| 
l + 9 v 



rank[G] 



(34) 



Notice that the first term in the r.h.s. is the central charge of a WZW model for the algebra 
G at level I (g v is the dual Coxeter number). It is now easy to check that the Y-system (|31I32|) 
becomes equivalent to the C^+i-related system of [17] at I = N, provided the infinite TBA 
diagram is truncated in the horizontal direction up to the node k (see figure and the following 
identifications are made 



Y°(9) = l/Y p a [ KN l (6) , g v = k + 2 , rank[G] = k + 1 , 



and hence 



7T 



\ i + ya J 

all nodes P 



iV(A; + l)(2(fc + !) + !) 
N + k + 2 



k — 1. 



(35) 
(36) 
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After the elimination of the massive nodes instead the correspondence is with a system of A k 
type at level I = 2N. The second contribution to (|33j) is 

<* = ? £ = 2w+¥+i ~ ' (37) 

magnonic nodes y 

Notice now that 

ci - c 2 = c[SU(2N) k+1 ) - c[S P (N) k+1 ] , (38) 

is equal to the difference in the central charges of the S , C/(2A^)fe +1 and Sp(N)k + i WZW mod- 
els. It is then natural to conjecture that the level k truncated system describes an integrable 
perturbation of the conformal model SU(2N)i c +i/ Sp(N)f : +i- The perturbing operator is easily 
identified by noticing that the solutions of the Y-systems fulfil the periodicity property [14, 24] 

Y^e + i,^^j=Y p ^(9), (39) 

then the argument of [14] leads to the conformal dimensions 

2N 

h = h = l- — , (40) 

N + k + 2 ' v ; 

suggesting that the models in question are most certainly perturbed by the operator (1, adj). 
On the other hand, in the k — > oo limit h = h — > 1 and 

c uv = lim (ci - c 2 ) = 27V 2 - iV - 1 , (41) 

k^oo 

which is the expected result (|25|) for the All nonlinear sigma model. 



6 Conclusions 

The main purpose of this paper is to show that the SU(2N) Gross-Neveu S-matrix, with a 
CDD factor selecting only even fundamental highest weight multiplets, correctly describes the 
ultraviolet limit of the sigma models on the All symmetric spaces. Considering also the earlier 
result of [12], where the (^2,1^2) scattering amplitude was checked by comparing the free energy 
in strong external field with perturbation theory, we can confidently confirm the correctness of 
the proposal [12]. 

Another result of the paper is the mapping, through the Y-system, between our models and 
particular limits of systems proposed by Kuniba and Nakanishi in [17]. This led us to the in- 
teresting conjecture that the truncated TBA equations describe the coset SU{2N)k+i/ Sp(N)k+i 
models perturbed by the operator (1, adj). As a follow-up project it would be interesting to check 
the integrability of these perturbed conformal field theories, to propose the exact S-matrix and 
to rigorously derive the set of thermodynamic equations. 

A more general conclusion is again related to the interesting fact that the finite-size physics of 
these more exotic models is again described by set of equations of the kinds [17,21]. It would be 
nice to understand this correspondence at a deeper level and to find the complete classification 
of the thermodynamic equations for integrable sigma models on symmetric spaces. 
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7 Appendix 

The definition of Fourier transform is 

—e^f(u) , /(a;) = J d6 e~^f(6) , (42) 

where g = 2N is the dual Coxeter number of SU(2N). The kernels A^\lo) and ^2a2b ar e 



fWc ^ o ut ^ sinh(a; min(a, 6)) sinh(u;(-L — max(a, 6))) 

sinh(wL) sinh(o;) 



^a6 M = 2cosh(a;) :nT7: A , (43) 



? _ M smh(2a;min(a,&))sinh(2a;(/V - max(a,6))) 

V 2a 2fe(w) = e ■ urn at\ ■ U7 \ ' i 44 i 

sinn (2a; TV) sinh(w) 

^>=(^V)r=^-^j), (45) 

a p (a>) = e~ pM . (46) 

j(L-i) _ 2(5 afe _ C^' 1 and C^ 1 - -1 are, respectively, the incidence and Cartan matrices of the 
algebra Al-i- The following relations are also useful 

CO <- 

6=1 V ; 

AT-1 



6=1 

AT-1 



^ ^ ] (2o;)m 26 cosh 6 = , (48) 

6=1 

E *2" m - ( 2 + C) ■ m 

In equation (|48|) m<ib = Msin(^) is the Perron-Frobenius eigenvector of the A^-i Cartan 
matrix. 
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